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Q ■ Abstract 

f^ , We describe how Computational Group Theory provides tools for manip- 

^^ I ulating tensors in explicit index notation. In special, we present an algorithm 

PLc that puts tensors with free indices obeying permutation symmetries into the 

(— I , canonical form. The method is based on algorithms for determining the canon- 

ical coset representative of a subgroup of the symmetric group. The complexity 
of our algorithm is polynomial on the number of indices and is useful for imple- 
mentating general purpose tensor packages on the computer algebra systems. 
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1 Introduction 

The connection between Group Theory and Tensor Calculus was established a long 
time ago. For instance, Weyl [|1| showed the relation of the symmetric group and 
group rings to tensor symmetries and tensor expressions respectively; Littlewood 
developed tools in group representation theory to address the problem of determining 
the dimension of the space generated by tensor monomials. Recently, Fulling et 
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al. used these tools to determine the number of independent monomials up to 
order 12 built out of the Riemann tensor and its covariant derivative. On the other 
hand, these authors have not established a method to determine the independent 
monomials explicitly. The present work performs an important step toward this 
direction. 

The application of Group Theory to develop algorithms for tensor manipulation 
is addressed in refs. [|| and 0. Ref. |^ uses group algebra which generates an 
algorithm of exponential complexity that can address tensor expressions with at 
most 11 or 12 indices. Practical applications demand better results. Ref. uses 
the backtrack algorithm to find the canonical form of tensor expressions built out of 
totally symmetric or antisymmetric tensors. 

The present work addresses the problem of finding efficient algorithms for ab- 
stract tensor manipulation using Computational Group Theory. The main problem 
consists in simplifying tensor expressions, which can be solved if one knows an effi- 
cient algorithm that puts tensor expressions into the canonical form. Ref. shows 
that the problem of finding the canonical form of a generic tensor expression reduces 
to finding the canonical forms of single tensors. At this point, the group-theoretic 
approach is the natural language to express the problem, since all informations about 
a single tensor can be represented by Group Theory, and can be efficiently processed 
using Computational Group Theory. 

We suppose that the tensors have only free indices and obey what we now call 
permutation symmetries, which are a set of tensor equations of the form 

'T'il ■■■ in ^ , n^ cr(ii ■■■in) /"i \ 

where a{ii ■ ■ ■ i„) is a permutation of zi ■ ■ ■ in and e^ is either 1 or —1. We present 
a polynomial-time algorithm to find the canonical forms of these tensors. This algo- 
rithm is a straightforward extension of algorithms to find the canonical coset repre- 
sentative of a subgroup of the symmetric group. 

In section 2 we describe the representation theory for tensors. In section 3 we 
present the main algorithm and discuss its complexity. 

2 Representation theory for tensors 

In this work we use the abstract-index notation for tensor expressions such as de- 
scribed in Penrose and Rindler's book [0]. We take Lovelock and Rund [§ as a 
general reference for Tensor Calculus. Ref. describes in details how a generic ten- 
sor expression can be converted into a sum of single tensors. To sum up, a generic 
tensor expression can be expanded so that it is a sum of tensor monomials. Each 
monomial is merged into a single tensor that inherits the symmetries of the original 
tensors. If each single tensor can be put into the canonical form, then the original 
tensor expression can be put also. At this point the problem consists in finding the 
canonical form of a single tensor obeying permutation symmetries. 



In the present context, we need only three kind of informations about a tensor: 
sign, index configuration, and symmetries. For example, if one wants to find the 
canonical form of T'^''"' knowing that the rank-3 tensor T is totally antisymmetric, 
one starts with 

{+l,[c,6,a]} 

and ends up with 

{-l,[a,b,c\}. 

The natural canonical configuration is —T°-'"^. 

In order to use a group-theoretical approach for this kind of manipulation we 
have to represent the symmetry as some group which acts on the index configuration 
and on the sign. This goal is achieved in the following way. 

Definition 1 (Symmetry of rank-n tensor from the group-theoretic point of view.) 
Let Sn be the symmetric group on the set of points {1, 2, ■ ■ ■ , n} and H the group 
({+1,-1}, X ) with multiplicative operation. A tensor symmetry S" is a proper 
subgroup of the external direct product H ® Sn such that (—1, id) ^ S, where id is 
the identity of Sn- 

In order to fix notation, from now on we assume that H and S are the groups 
described in Def. 1. We also assume that G is the subgroup of Sn such that H®G =< 
(—1, id), S >, i.e. H ® G is the smallest subgroup oi H ^ Sn which contains S and 
{—l,id). Therefore, the order oi H ^ G is two times the order of S. 

Now let us describe formally how S acts on a tensor. Each element of 5* is a pair 
consisting of a sign (±1) and a permutation. The sign of a permutation n will be 
denoted by e^. So, an element of S has the form (e^, tt), where e^, = ±1 and n E G. 
The action of s on a totally contravariant rank-n tensor with index configuration 



frpili2---in\ __ , rpiiTT i27r ■■■i„7r /q') 

where the notation ii^ means that the subscript of i is the image of point 1 under 
the action of permutation tt. This notation clearly shows that the permutation acts 
on the positions of the indices and seems to be superior to the one used in (|T]). 

The notation +7r stands for (1, tt) and —n for (—1, vr). For example, the group 5" 
which describes the symmetry of a totally antisymmetric rank-3 tensor is 

S={+td, -(1,2), -(1,3), -(2,3), +(1,2,3), +(1,3,2)}, (3) 

which, in tensor notation, corresponds to 

rpabc rpabc rpabc rpbac rpahc rpcba 

rpabc rpacb rpabc rpbca rpabc rpcab / a\ 



We use the term "permutation" to describe the elements of S extending the syntax 
used for elements of S'„. The term "permutation sign" refers to e^r- The sign e^ has 
no relation to the parity of it in the general case. 
With no loss of generality, we take 

as the standard configuration. All other configurations are obtained by acting per- 
mutations oi H ^ Sn on (^. A generic configuration is denoted by T^^'"^". We use 
the sequence of j's (ji ■ ■ ■ j„) as a generic permutation oiii- ■ ■ in- The standard con- 
figuration (^ is associated with the element {+l,id), which is the minimal element 
oi H ® Sn- 

In Def. 1, the requirement that {—l,id) ^ S avoids the cancelation of a tensor 
independently of its components. Note that if both (+l,7r) and (— l,7r) are in 5*, 
then we come to two conclusions: there are two equal index configurations such that 
ThJ2-i^ = _Tiii2-jn and {-l,id} G S. Vice-versa, if {-l,id) e S then (-1,^) G S 
if and only if (+1, tt) G S. When implementing a tensor package in some computer 
language, it is useful to have an efficient method to exhibit the symmetries that can- 
cels the tensor independently of its components. Algorithms for testing membership 
provide such method. 

The index configuration T^^^^'"^" is equivalent to T^^^^'"^" if and only if there is 
an element in S such that 

In other words, there is an element (+l,7r) G S such that the lists J1J2 ■ ■ ■ jn and 
iiTT 12^ ■ ■ ■ in^ are exactly the same. 

Proposition 1 The set of index configurations equivalent to T*^*^'"*" is given by the 
action of S over T*^*^"*". The cardinality of this set is the order of G (|G|). 

Note that \G\ = \S\. 

Proposition 2 Let j">-i^'t2^---^n^ Iqq an index configuration which is not equivalent to 
±T*i*2---«i (^i.e. both +7r) and — tt) are not in S). The set of configurations equivalent 
^o T*!" *2T---«„7r jg ^]-^g right coset oi S in H ® Sn which contains (+1, tt). The number 
of index configurations equivalent to J'«it«2t---«„t jg j^^j^ 

The final part of proposition 2 follows from the fact that all cosets oi S in H ® Sn 
have cardinality |G|. 

Proposition 3 Consider a set of independent index configurations with positive 
coefficients. The maximum cardinality of this set is the index of G in Sn {\Sn'-G\). 



Usually, no one describes the symmetry of a tensor in index notation by listing 
all equivalent index configurations. In general, one gives a few equations and skips 
the ones that can be obtained from the original equations. The following example 
illustrates this situation. Let T"'"^'^ be a rank-4 tensor with the following permutation 
symmetries. 

rpabcd rpbacd ('7a\ 

rpabcd rpcdab ('7\\\ 

iFiom. these equations we can obtain 

rpabcd rpabdc fQ\ 

Note that one can describe the symmetries of T"'"^'^ by using a different set of equa- 
tions, for example eqs. (|7|) and (|^). In group-theoretic language, this is equivalent 
to describe a group by generators. In the example above, the generating set K for 
the symmetry of T"'"^'^, described by eqs. ([^), is 

K = {-(1,2),+(1,3)(2,4)}. (9) 

The symmetry S is the group generated by K {<K>) 

S = {+id, -(1,2), -(3,4), +(1,2)(3,4), +(1,3)(2,4), 

-(1,32,4), -(1,4,2,3), +(1,4)(2,3)}. (10) 



The group G is generated by {(1, 2), (1, 3)(2, 4)} and is obtained from (|T0|) by remov- 
ing the permutation signs. /^From proposition 1, the set of configurations equivalent 
to the standard configuration T"'^'^'^ is given by the action of all elements of S on 
rpabcd^ which yields 

f rpabcd rpbacd rpabdc rpbadc rpcdab rpcdba rpdcab rpdcba\ ('\'\\ 

The number of equivalent configurations is 8, the order of 5". Now consider T"'^'"^, 
which is not in set (|11]). This index configuration is obtained from T"'^'^'^ by the action 
of +(2,3). Neither +(2,3) nor —(2,3) are in S. Proposition 2 says that the set of 
configurations equivalent to T"'^''"' is the right coset oi S m. H ® S^ which contains 
+ (2, 3). This coset is obtained by multiplying each element of S by +(2, 3): 

S X (+(2,3)) = {+(2,3), -(1,3,2), -(2,3,4), +(1,3,4,2), 

+ (1,2,4,3), -(1,2,4), -(1,4,3), +(1,4)}. (12) 

The action of this coset on T"'^'^'^ generates all index configurations equivalent to 
rpacbd^ Note that |S'4:G'| = 3. From proposition 3 we know that the set of independent 
configurations has cardinality 3. An example of this set is 

J rpabcd rpacdb rpadcb\ ('\'X\ 



which is a complete right transversal for G in S^ in tensor notation. One can 
easily recognize that eqs. (|^) and (|^) represent the symmetries of the Riemann 
tensor without taking into account the cyclic symmetry. The group generated by 
{(1,2), (1,3)(2,4)} is the dihedral group of order 8 {D^). Then, the symmetry of 
the Riemann tensor given by ([1^) is the largest subgroup oi H ® D^ which does not 
contain (— 1, zd). 

In order to address the problem of finding the canonical configuration equivalent 
to a given index configuration, we have to define an order for the permutations of 
Sn- Let b = [6i, ■ ■ ■ , 6„] be a list of n distinct points of the set {1, ■ ■ ■ , n}. Suppose 
that pk and pi are points. We define the order "-<5" for points with respect to b in 
the following way: pk -< Pi H k < I. So, p^ is smaller that pi if p^ comes before pi in 
b. We omit the index b from "^b" to simplify the notation, b^ = [bl,- ■ ■ , 6^] is the 
image of b under tt. Define 

C = {b\ TV e Sn}. (14) 

Now let extend the order "-<" to the elements of C. Let Li and L2 be in C. Li -< L2 
if 

Li[l]^L2[l] or (Li[l]=L2[l] and Li[2..n] ^ L2[2..ri]), (15) 

where L[i] means the i-th element of L and L[2..?7,] means {L[2], L[3], ■ ■ ■ , L[?t,]}. If 
TTi and 7r2 are in Sn then 

TTi ^ 7r2 ^ 6^^ ^ 6^^ (16) 

Now we extend the order "^" to group S: 

(e^,,7ri) -< (6^2,712) -^ TTi -< 7r2. (17) 

Recall that (+l,7r) and (— l,7r) cannot be at the same time in S. So we simply 
disregard the sign as stated in (|T7|) . "^" is a well-order (total order with a minimal 
element) in Sn, S, and in any coset of S* in if ® Sn- /,From now on, "minimal point" 
refers to the order "-<" with respect to some 6. 

A canonical right transversal for G in Sn {S in H ^ Sn) is a complete right 
transversal for G in Sn {S in H ® Sn) such that each coset representative is minimal. 
The set ([T3|) is the canonical right transversal in tensor notation for group G in 6*4 
with respect to b = [1,3,2,4], where G is generated by {(1,2), (1,3)(2,4)}. Note 
that the order "-<" allows to sort a canonical right transversal for G in Sn- This is 
important for addressing the simplification problem when there are side relations, 
such as the cyclic symmetry. This kind of symmetry is not addressed in this paper. 

3 Algorithm to canonicalize tensors with free in- 
dices 

Now we address the following problem. Suppose one gives a totally contravariant 
tensor with the symmetries described by a set of tensor equations, and a free index 
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configuration that is not the standard one. Find the canonical index configuration 
with respect to the order "-<" . For example, suppose that rank-4 tensor T has the 
symmetries (|^), and one gives the following index configuration: T''™'^. What is the 
canonical configuration with respect to h = [1,3,2,4]? 

Using the representation theory, the problem above can be solved if one knows 
the solution of the following problem. Given a generating set K for the group S 
and an element (e^r, tt) in i7 ® Sn, find the canonical coset representative of the coset 
S X (e^, tt) with respect to the order defined by (\i7\} . 

The answer to this problem is an algorithm that can be used to put a tensor with 
free indices into the canonical form. Before describing the algorithm, we extend some 
well known definitions for permutations groups, in order to use them in connection 
with S, which is a direct product of groups. The extensions are straightforward. 

Definition 2 Let p be a point in the set {1, ■ ■ ■ ,n}. The stabilizer of p in S is the 
subgroup Sp defined by 

Sp = {seS\f = p}, (18) 

where p^ = p'^ and s = (e,r, tt). 

In other words, Sp consists of all elements of S that fix the point p. Sp is a subgroup 
of S. If s = (e^,7r) G S, "s fixes p" means "tt fixes p" . 

Definition 3 Let Q be a subset of the set of points {1, ■ ■ ■ ,n}. The pointwise 
stabilizer of Q in S is the subgroup Sq defined by 

SQ = {seS\WqeQ,q' = q}. (19) 

In other words, Sq fixes all points of the subset Q. 

Definition 4 A ordered subset b = [6i, ■ ■ ■ , bm] of the set of points {1, ■ ■ ■ ,n} is a 
base for S if Sbj^^...^b,n = (l; ^^)- 

This means that the only element of S that fixes all points of b is the identity. A 
useful property of a base is that an element of S is uniquely determined by the base 
image. We can order the points {1, ■ ■ ■ ,n} such that the base points are the first 
m points. Let us name the remaining points /i, ■ ■ ■ ,ln-m- They simply follow the 
usual increasing order. From now on, the order "-<" defined in the previous section 
is based on the set [6i, ■ ■ ■ , bm, h,- ■ ■ , In-m]- 

Definition 5 A strong generating set K for S relative to the base 6 = [6i, ■ ■ ■ , bm\ 
is a generating set for S with the following property: K fl Sh^^... ^hj is a generating set 
for the pointwise stabilizer Sbi,---,b ■, for 1 < j < m — 1. 



In other words, a generating set K for S is strong if, after selecting the permutations 
of K that fix the points 61, ■ ■ ■ , hj, one has a set that generates the group S'b^...._b.. 
This must be vahd for j from 1 to m — 1. 

The input of the algorithm is 

(a) an index configuration T-'i-'^ ■"■'"; and 

(b) a base and a strong generating set for S. 

The output is the canonical index configuration, which is obtained by the action 
of the canonical coset representative on the standard configuration r*^*^"*". If S is 
described by a generating set that is not strong, the first step is to obtain a strong 
generating set and a corresponding base. Refs. and [1^ present algorithms that 
perform this task. These algorithms must be extended in order to work within 

H ®Sn. 

The algorithm Canonical described below uses the general structure of chain 
of stabilizers developed by Sims |]TT| and is an straightforward modification of the 



algorithm presented by Butler |^|. Following Sim's notation, let 5''-*^ be the group 
that stabilizes the points {6i,--- ,6j_i}, i.e. S"*^*^ = Si,^^... ^b^_^. Then S^^^ = S and 
^M = (1, id). Let A(*) be the orbit of S^'^ that contains the point k, i.e. A^*) = bf\ 
Suppose that i^ is a strong generating set for S and define K^'^'> = K H S^'^K K^''^ 
is a strong generating set for S^^\ for 1 < i < m. Let u^^^ be the Schreier vector of 
A^*) with respect to generators K^'^K Here, the components of the Schreier vectors 
are elements of S. Following Butler [|10|, if g G A*^*\ let trace{q, u^^^) be the element 
{euj,io) of S^^^ such that p'^ = q, where p is the minimal point of the subsets of A*^*^ 
that contains q. 

The algorithm consists of m loops. Suppose that [gi, ■ ■ ■ , Qm] is the base image of 
the canonical representative. The i-th loop finds permutation (eA, A) that determines 
qi, i.e. bf = qt. The permutation (eA, A) obeys the constraint 

[(6i)\---,(6._i)^] = [gi,---,g._i]. (20) 

The set of all elements of S that obey (^) is given by S^^^ x (eA, A), then b^ ^ 
yields all possible images of bi in the coset S^^^ x (ca, A), qi is the minimal point of 
these images. At the last loop, (ca, A) gives the complete base image [qi,- ■ ■ , qm]- 



Algorithm Canonical (free indices) 

Input: T^'i-J" = g^7^n-i2— in-^ ^here (e^,7r) G S; 
b = [bi,- ■ ■ , bm] base for S; and 
Ks strong generating set for S with respect to b. 



Output: caT*!^ *2^"*n^j where (e^, A) is the canonical representative of the 
coset of S in H ^ Sn that contains the permutation (e^, tt). 

begin 

(* initiahzation of (eA, A) and K *) 
(eA,A) := (e^,7r); 
K := Ks; 

for i from 1 to m do 

(* A is the basic orbit A*^*) *) 

A := 6<^>; 

k := position of the minimal point of A'^; 

p := k-th point of A; 

{e^^jUj) := trace{p, u), where u is the Schreier vector of A with respect to K; 

(eA,A) := (e^.u) x (eA,A); 

K := remove permutations of K that have point bf, 
end for; 

return eAT*i^*2^'"*"^; 
end 



Let us see an example. If one gives the index configuration T''™'^, where T has 
the symmetries (|^), the element {e'^,Tc) is +(1,2,3). A base for S is [1,3], so the 
order "-<" is based on the list b = [1, 3, 2, 4]. A strong generating set with respect to 
base [1, 3] is 

K = {-(l,2),-(3,4), +(1,3)(2,4)}. (21) 

First loop yields: A := {1,2,3,4}; A^^'^'^) ;= {2,3,1,4}; k := 3; p := 3; (e^,cj) : = 
+ (1,3)(2,4), since 6^ = p; (ca, A) := +(2,4,3); K := {-(3,4)}. (eA,A) applied on 
rpabcd gj^gg r^adbc_ Second loop yields: A := {3,4}; A^^'^.s) .= {2,3}; k := 2, since 
3 ^ 2; p := 4; (e^,a;) := -(3,4); (eA,A) := -(2,4); K := {}. The algorithm finishes 
and the canonical configuration is eA T*i^*2^*3^*4^ = —j'adcb^ 

The algorithms to find strong generating set, basic orbit, Schreier vector, and 
trace for permutation groups are described in refs. |^ and ||10|- The extension of 



these algorithms to work within the direct product H ^ Sn is straightforward if one 
uses the fact that a image of a point p under the action of (e^r, vr) is p'^. The product 
of permutations in Sn is naturally extended to the product of elements of i/ ® Sn- 

The analysis of the complexity of the algorithm Canonical is the following. Al- 
gorithms to find basic orbit, Schreier vector and trace have an O(n^) bound. Since 
the algorithm Canonical performs n loops in the worst case, it is bounded by 0{n^). 
A strong generating set for symmetry S is required. It is known that Schreier-Sims 
algorithm has an 0{n^) bound. So, if the generating set of S is not strong, the overall 
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bound for the algorithm to find the canonical form of tensors is 0{n^), where n is 
the number of indices. 
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